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We propose two methods of estimating a systematic error in extrapolation to the infinite-size 
limit in the study of measuring the Haldane gaps of the one-dimensional Heisenberg anti- 
ferromagnet with the integer spin up to S = 5. The finite-size gaps obtained by numerical 
diagonalizations based on Lanczos algorithm are presented for sizes that have not previously 
been reported. The changes of boundary conditions are also examined. We successfully demon- 
strate that our methods of extrapolation work well. The Haldane gap for S = 1 is estimated to 
be 0.4104789 ± 0.0000013. We successfully obtain the gaps up to 5 = 5, which make us confirm 
the asymptotic formula of the Haldane gap in S oo. 
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1. Introduction 

Extrapolation is a fundamental technique in a lot of 
studies in physics. In particular, the technique is of- 
ten carried out in the condensed-matter physics when 
one attempts to know a quantity in the thermodynamic 
limit from several finite-size data. One of the reliable 
ways to obtain such finite-size data is the numerical- 
diagonalization method applied to the Hamiltonian ma- 
trix describing a system. This method provides us with 
very precise finite-size data although available system 
sizes are limited to being very small. This method is 
non-biased against the effects of interaction; thus it con- 
tributes much to the understanding of many-body prob- 
lems. Typical examples are the quantum spin systems, 
in which there often appear nontrivial quantum states 
due to the presence of interactions between spins. One 
of them is the ground state of the integer-spin one- 
dimensional Heisenberg antiferromagnet. In this system, 
an energy gap exists between the unique ground state 
and the first excited state; this gap is called the Haldane 

1 2 

gap.^'^ 

The magnitude of the Haldane gap has been estimated 
by various numerical methods."^ There are three represen- 
tative approaches. The first one is the numerical Lanc- 
zos diagonalization of finite-size clusters. In the S — 1 
case, system sizes up to 22 sites were treated under the 
periodic boundary condition."' Since the available sys- 
tem sizes are small, an appropriate extrapolation is re- 
quired. Unfortunately, it is difficult to estimate a system- 
atic error due to the extrapolation. The second one is the 
density matrix renormalization group (DMRG) method.^ 
The calculation was carried out under a peculiar bound- 
ary condition, namely, each edge of the S = 1 chain con- 
necting with an S* = 1/2 spin. In this way, it is necessary 
to tune the artificial interaction at the edges. The third 
one is a quantum Monte Carlo (QMC) simulation.^ Since 
this simulation was performed by the loop algorithm to- 
gether with a continuous imaginary time technique, cal- 
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culations of very large systems are available. However, a 
statistical error due to a Monte Carlo sampling cannot 
be avoided. Up to the present time, these approaches give 
consistent estimates of the Haldane gap with their own 
errors of the same order. 

Under such circumstances, we attempt again to es- 
timate the Haldane gaps of the Heisenberg antiferro- 
magnetic spin chain as precisely as possible. Our main 
method is the numerical diagonalization. A primary pur- 
pose of this paper is to propose a procedure to obtain a 
reliable error in extrapolation of finite-size data toward 
the infinite size of the system. 

In the extrapolation, the weak system-size dependence 
of the finite-size gap is favored. The system-size depen- 
dence of the gap is determined by the choice of the 
boundary condition. It is known that the dependence 
becomes suppressed when one twists the boundary con- 
dition from the periodic boundary condition.^ Under this 
background, we examine boundary conditions in the case 
of 5' = 1 to know which boundary condition is appropri- 
ate. We then find that the twisted boundary condition is 
the most appropriate one among periodic, twisted, and 
open boundary conditions. The twisted boundary condi- 
tion gives a good sequence of finite-size excitation gaps 
for the sake of the extrapolation. Next, we develop a 
method to obtain a reliable error when the convergence of 
the data sequence can been accelerated. Thereby, a very 
precise estimation of the S = 1 Haldane gap is success- 
fully obtained. When one imposes the twisted boundary 
condition for 5' = 2, 3, 4, and 5, the Haldane gaps can 
be obtained to be nonzero values in the thermodynamic 
limit in spite of the fact that the Lanczos method can 
treat only the extremely small system sizes. When S be- 
comes larger, the convergence acceleration becomes more 
difficult. In such a case of the acceleration in failure, we 
also develop another procedure to estimate an error to 
the excitation gap of the infinite system. 

This paper is organized as follows. In the next section, 
the model Hamiltonian and the method of calculation 
will be explained. In the first half of §3, the numerical 
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results of the S = 1 system are presented and discussed. 
Boundary conditions are examined and the extrapola- 
tion based on the convergence acceleration is performed. 
We propose its error estimation and demonstrate the va- 
lidity. In the second half of §3, the cases of S > 2 are 
studied. Another procedure of obtaining an error is in- 
troduced. The final section is devoted to the summary 
and some remarks. 

2. Hamiltonian and Method 

The Hamiltonian of the present model is given by 

AT 

'm+1, (1) 



H = ^ JmS 



1=1 

where Sm is a spin operator with its amplitude S at site 
m. Here N is the number of spin sites. The system size 
N is supposed to bo an oven integer. The amplitude of 
the exchange interaction is denoted by Jm which will be 
defined later when the boundary conditions is explained. 

In this paper, we carry out numerical diagonalizations 
for finite-size clusters of systems of S' = 1, 2, 3, 4, and 5. 
We calculate the ground-state energy Eq and excitation 
energies (the first excitation Ei and the second excita- 
tion E2) by the method of Lanczos algorithm.* We have 
successfully developed a code for parallel calculations of 
the Lanczos algorithm. The maximum sizes in this paper 
are iV = 24 for 5 = 1, iV = 16 for 5 = 2, TV = 12 for 
5 = 3, TV = 12 for 5 = 4, and JV = 10 for S = 5. These 
sizes have not been treated in the Lanczos calculations 
as long as the present authors know. Note that we only 
assume the conservation of the ^-component of the total 
spin. Thus, arbitrary shapes of clusters can be treated. 
The dimensions of the Hilbert space are very large. For 
example, the dimension of the largest subspace, namely 
Sfot = 0, of TV = 24 for S' = 1 is 27 948 336 381, where 
S'tot is the z component of the total spin. It is worth to 
emphasize that the parallelization makes it possible to 
carry out the Lanczos calculations. 

3. Result and Discussion 

3.1 Case for S=l 

3.1.1 Boundary Conditions 

Let us consider differences from the choice of boundary 
conditions. The differences from boundary conditions ap- 
pear in results of calculations of finite-size systems. On 
the other hand, the differences are supposed to disap- 
pear in the limit of iV ^ 00. Finite-size effects depend 
on the choice of a boundary condition. Namely, a differ- 
ent type of boundary condition gives a different finite-size 
sequence concerning with a physical quantity. In order to 
obtain precisely the information in the thermodynamic 
limit that is not affected by boundary conditions, one 
should employ an appropriate boundary condition. Such 
a condition is not necessarily the periodic boundary con- 
dition. The appropriate condition depends on systems 
and physical quantities. Therefore, the examination of 
various boundary conditions is important. Here we focus 
our attention on the problem of the Haldane gaps and 
begin with such an examination in the 5=1 case. 
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Fig. 1. Energy differences of finite-size 5=1 systems. Open and 

closed circles denote results of Ei — Eq under the periodic bound- 
ary condition and the twisted boundary condition, respectively. 
Triangles and diamonds denote results of E2 — Ei and Ei — Eq 
under the open boundary condition, respectively. Note that the 
maximum system size is iV = 24 irrespective of boundary condi- 
tions. 



In this paper, we examine three types of the bound- 
ary conditions: the open, periodic, and twisted boundary 
conditions. The open boundary condition is given by 

(m < A'') 
(m = N). 

The periodic boundary condition is given by S^+i = 
Si and = 1 for arbitrary m. The twisted boundary 

condition is given by 



Jrr 



(2) 



N+1 



N+1 



-Sf, S'tv+i — Si , (3) 

and Jto = 1 for arbitrary m. Note here that energies are 
measured in units of nonzero Jm\ therefore we take it 
unity. 

First, we show numerical results of system size depen- 
dence of energy differences in Fig. 1. Our results under 
the periodic boundary condition agree with those in ref. 4 
up to A = 22; our results for A = 24 are new. One can 
clearly observe that E\ — Eq under the open boundary 
condition vanishes in the limit of A ^ oc'. This behav- 
ior is consistent with a quasi-degeneracy of the Haldane- 
type ground states. Under the open boundary condition, 
the Haldane gap appears above these degenerate ground 
states; the energy difference E2 — Ei decreases gradu- 
ally when A is increased and seems to converge around 
0.4. Under the twisted boundary condition, on the other 
hand, the energy difference Ei — Eq increases with in- 
creasing A. This dependence will also be very useful 
when we study gaps for the cases of S* > 2. In this sec- 
tion, let us compare the speed of the convergence of the 
finite-size sequence under each boundary condition. In 
order to achieve it, we consider the ratio defined as 

G«(A) - G«(A + to) 



i?«^(A,m) = 



Gi{N) - G^N + m)' 



(4) 



Here G''(A) is the energy difference with system size A 
under the C boundary condition which converges to the 
Haldane gap in the limit of A — > 00; namely, G^{N) = 
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Fig. 2. Ratio (Af, m) v.s. inversed system sizes. Circles denote 
the case of { = periodic, = open. Diamonds denote the case of 
{ = twisted, ^ = periodic. Closed symbols mean m = 10 which 
is the largest m in this figure. 



El — Eq for the periodic and twisted boundary condi- 
tions, and G'^=°P°"(A^) = E2 - Ei. Results of the ratio 
R^''{N, m) v.s. are depicted for two types in Fig. 2. 
All the absolute values of the presented ratios are less 
than unity. The ratio of ^ = periodic and C = open gets 
gradually smaller when N is increased. The ratio seems 
to vanish in the limit of iV ^ 00. If the ratio in the limit 
of m — > 00 vanishes and when each G^{N) converges to 
the same value, namely the Haldane gap A(S' = 1) in 
this case, one obtains 



lim 



G«(A^)- A(S'== 1) 



= 0. 



(5) 



G«(A^)-A(5 = 1, 

This suggests that the sequence under the periodic 
boundary condition converges faster than that under the 
open botmdary condition. (See appendix.) Thus the pe- 
riodic boundary condition is more appropriate than the 
open boundary condition to measure the Haldane gap 
within the system sizes that arc available in numerical 
diagonalization calculations. On the other hand, the ra- 
tio of ^ = twisted and C, = periodic in the limit of N ^ 00 
seems not to vanish but to converge to a nonzero value. 
This means that concerning with the convergence of the 
sequence, the speeds of cases of the twisted and peri- 
odic boundary conditions arc almost the same with each 
other. Even though the speeds are comparable, the ab- 
solute values of the ratios are much smaller than unity 
{\R^''{N,m)\ < 1/4). This suggests that each datum of 
the sequence under the twisted boundary condition is 
closer to the Haldane gap A(5 = 1) than the correspond- 
ing datum under the periodic boimdary condition. In 
this meaning, we can conclude that the twisted boundary 
condition is more appropriate than the periodic bound- 
ary condition. Therefore, the twisted boundary condition 
is the most appropriate among the present three condi- 
tions. 

3.1.2 Extrapolation 
In this subsection, we attempt to extrapolate the above 

sequences of the S = 1 finite-size energy differences by 
means of the technique of convergence acceleration. We 



here apply Wynn's epsilon algorithm^ given by 

1 1 1 

+ 



-(fc+i) _ Ak) 



A{k) _ Ak) ' .(k) _ .(k) 



l(fe) 



a 



Ak-l) _ Ak) ' 



(6) 



when a = 1. Here the initial condition is given by 
A^^^ = G{N) and A^^^^ = 00. It was in ref. 4 that 
this transformation was applied for the first time to es- 
timate the S = 1 Haldane gap. Note that when we take 
a = in eq. (6), the transformation is reduced to the the 
one called as the Aitken-Shanks process. The Aitken- 
Shanks transformation was used in ref. 11. Both of the 
transformations make us possible to accelerate the con- 
vergence of a finite sequence and to give a candidate for 
the extrapolated value. It is, however, difficult to obtain 
a systematic error only within the framework of each 
transformation. Under such circumstances, the authors 
of ref. 4 considered the variance of a within a successful 
acceleration of the convergence. We have found problems 
in the argument of estimating their systematic error in 
ref. 4 when we examine the convergence of the S = 1 
finite-size gaps up to A'' = 24 under the periodic bound- 
ary condition. We present the results of the table of the 
convergence for a = 1 in Table I. Note here that a part 
up to A' = 22 in Table I was reported and that the en- 
ergy difference of A = 24 and its posterity are new. Let 
us mention the problems while we are reviewing the pro- 
cedure of ref. 4. 

Let us explain the procedure of ref. 4 briefly. First, one 
considers the following decay lengths defined as 



4^) =2/ log f^-- 



A 



(k) ^ 

N-2 



A 



(fc) 



A 



(k) 



(7) 

'■N-2 ^'N ) 

in order to examine converging behavior of the sequence 
for each step k. One should note that increases 
monotonically. This means that the convergence becomes 
slower when A^ is getting larger. This is a source of dif- 
ficulties in estimating the extrapolated gap. In order to 
overcome this difficulty, an acceleration is introduced. To 
examine whether the acceleration of yl]^ for each k from 



A]y is successful or not, ref. 4 investigates the follow- 
ing three conditions. 

I a!>^ for each k is monotonic. 

II for each k is monotonically increasing, namely, 
the following condition holds, 

SJV-l-2 ^ ?JV • 

III The following condition holds, 

SAT ^ SJV • 



(8) 



(9) 



Conditions I and II suggest that properties of the initial 
sequence are preserved even though the acceleration is 
carried out. Condition III means whether an element 
with a long decay length is successfully removed by the 



acceleration, ref. 4 considered that A 



(5) 

N=22 



0.410498 



is successfully accelerated and that it is reliable as an 
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Table I. Sequence of finite-size gaps for the 5=1 case under the periodic boundary condition and the convergence acceleration based 
on Wynn's algorithm (cq.(6) with a = 1). 



N 


r;(M\ — 4(0' 


AO) 




Ai) 




A^) 


.(3) 


A-i) 






.(5) 


2 

4 

a 
u 

8 


2 

0.593555254375 


1 ^7 
2.54 


n fil 9*^909'^ 

0.48753251 


















10 


0.524807950414 


3.26 


0.44377567 


1.91 


0.43525877 














12 


0.484196469912 


3.80 


0.42557752 


2.28 


0.41798489 














14 


0.458965346938 


4.20 


0.41757427 


2.43 


0.41308940 


1.59 


0.41258393 










16 


0.442795561359 


4.50 


0.41394088 


2.53 


0.41152416 


1.75 


0.41114647 










18 


0.432221469865 


4.71 


0.41223978 


2.64 


0.41095437 


1.98 


0.41074380 


1.57 


0.41071215 






20 


0.425210314459 


4.87 


0.41141374 


2.77 


0.41071416 


2.32 


0.41058977 


2.08 


0.41055478 






22 


0.420515020390 


4.99 


0.41099554 


2.94 


0.41059985 


2.69 


0.41052336 


2.38 


0.41050133 


1.85 


0.41049811 


24 


0.417346883838 


5.08 


0.41077448 


3.14 


0.41054158 


2.97 


0.41049612 


2.24 


0.41048618 


1.59 


0.41048426 



approximate for the gap value. In ref. 4, the region of a 
where all the three conditions hold was found around a = 
1 within data up to A'^ = 22 even when one applies the 
acceleration iteratively up to /c = 5. Finally, the obtained 
region of a gave a systematic error. 

Let us examine the situations in all the presently avail- 
able data up to = 24. In Table I, up to /e = 2 
are monotonically increasing. However, one finds that 

(3) 

does not show a monotonic A''-dependence and that 

S}^^ decreases with N . Condition (8) docs not hold. This 



new sequence defined as 



means that it is unclear whether the acceleration of A 



(3) 



, and A^^ are successful or not. Next, we examine 
the behavior of and when we tune a around 

f3) 

a = 1. The decreasing behavior of disappears around 
a. = 1.2; however, s}^^ is still decreasing. The present ex- 
amination of data up to iV = 24 suggests that it becomes 
unclear whether ^^^^22 — 0.410498 is appropriate as a 
reliable estimate or not according to the criteria of the 
above three conditions. At least for a = 1, the accelera- 
tion of A^^^ up to fc = 2 seems to be successful. It is pos- 
sible to use A^^ instead of ^^^22 ^ reliable estimate. 

(2) 

In this case, let us remember that A)^ is monotonically 

(2) 

decreasing. This suggests that any of data in A]^^ gives 
an upper bound for the true Haldane gap. Even though 
the value of q is tuned, one obtains only an assembly of 
upper bounds. There is no evidence to show that the true 
gap value is in the region of the assembly. Thus, tuning 
a is not an appropriate way to obtain a reliable error of 
the Haldane gap. Therefore, we have to develop another 
strategy without tuning a to achieve it. 

Next, we present the result of the convergence acceler- 
ation of our gap data under the twisted boundary condi- 
tion. The table for a = 1 is shown in Table II. One can 
observe that all the above conditions hold. It is reason- 
able to consider that all of A'^' are successfully accel- 
erated. The difference between the periodic and twisted 

(k^ (k) 

boundary conditions is the direction of A)^^'; namely A)^' 
in Table II is monotonically increasing. This means that 
the data in A^jl^^ are lower bounds. In order to obtain a re- 
liable error only from the data under the twisted bound- 
ary condition, it is required to create another sequence 
that is monotonically decreasing. Here let us consider a 



^"""^ 'S'jv-i-2 — Sn — T]sr+2 + Tj 



N 



(10) 



where Tjv = and Sn = A^J^ If A^JP is success- 



fully accelerated from , it is expected that the se- 

(k) 

quence B]v+i convergent from the side opposite to 



T/v = A^^'^ and Sn 



A^j^ ^\ (See appendix.) In Ta- 



ble III, we present the result of obtained from Aj^' 
under the twisted boundary condition. One observes that 
all of B^-* are decreasing with increasing N. The differ- 
ence of the dependences of and B^j^^ makes us know 
that the true gap value is between them. In this works, we 
employ A^-* and B^-* whose dependences are confirmed 
to be opposite to play safe. Namely, the 5 = 1 Hal- 
dane gap is expected to be between ^^^^24 = 0.41047777 
and B^^}^23 — 0.41048023. Therefore, our present conse- 
quence for the Haldane gap for 5 = 1 is 

0.4104789 ±0.0000013. (11) 

This estimate agrees with the estimate from the QMC 
method^ and that from the DMRG one.^ Note that our 
estimate (11) is more precise than any other estimates 
as long as the present authors know. Before finishing 
this paragraph, we illustrate and b'^^ in Fig. 3 so 
that one visually captures the features of these sequences, 
which are explained in the above. 

The sequence is a general way to estimate a re- 
liable error; the case of the periodic boundary condition 
in Table I is applicable. In this case, is created from 

Sn = Ajv'' ^'^'^ 



A^^'' because it is confirmed that 
A^^^ is successfully accelerated from A^j^\ Note that in 
this case, B^j^^ is monotonically increasing and converg- 

(k) 

ing from the smaller side. Thus, B)^' gives a lower bound 
of the gap value. The data up to = 24 under the peri- 
odic boundary condition suggests that the Haldane gap 



is between ^^^23 and ^^^24! ^e have 0.41050 ±0.00005 
as an estimate in the periodic boundary condition. This 
result is also consistent with the estimate (11), the one 
from the QMC method,^ and the one from the DMRG 
one.^ This indicates that the present procedure makes it 
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Table II. Sequence of finite-size gaps for the 5 = 1 case under the twisted boundary condition and the convergence acceleration based 
on Wynn's algorithm (cq.(6) with a = 1). 



N 


Lz(l\ ) — /Ijy 


AO) 




Ai) 




A^) 


.(3) 


A-i) 






.(5) 


4 


n 9Q77fiQ^7Q1 ^1 






















a 
u 

8 


u.ouzu J-O^yoo J-O 
0.386237672978 


1.98 


0.39977728 


















10 


0.396943190982 


2.53 


0.40581471 


















12 


0.402443823536 


3.00 


0.40825701 


2.21 


0.40920572 














14 


0.405509288158 


3.42 


0.40936819 


2.54 


0.40998448 














16 


0.407315632794 


3.78 


0.40990703 


2.76 


0.41027702 


2.04 


0.41035307 










18 


0.408423139414 


4.09 


0.41017830 


2.91 


0.41039295 


2.16 


0.41043554 










20 


0.409122144092 


4.35 


0.41031824 


3.02 


0.41044086 


2.26 


0.41046285 


1.81 


0.41046801 






22 


0.409572951736 


4.56 


0.41039177 


3.11 


0.41046146 


2.37 


0.41047263 


1.95 


0.41047535 






24 


0.409868488828 


4.74 


0.41043100 


3.18 


0.41047068 


2.49 


0.41047640 


2.10 


0.41047777 


1.80 


0.41047815 



(fc) 

Table III. Upper-bound sequence B)^' of the 5 = 1 Haldane gap. 



N 










9 


0.41728860 








11 


0.41289925 








13 


0.41156232 


0.41142826 






15 


0.41100869 


0.41071649 






17 


0.41074765 


0.41055312 


0.41054058 




19 


0.41061764 


0.41050468 


0.41049203 




21 


0.41055134 


0.41048859 


0.41048271 


0.41048358 


23 
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0.41048286 


0.41048036 


0.41048023 
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Fig. 3. The system size dependence of the sequences and 
-B]y under the twisted boundary condition. Open circles, tri- 
angles, squares, pluses denote A^^\ A^^\ and A^^\ re- 
spectively. Closed circles, reversed triangles, diamonds, crosses 
denote B^j^\ ^n^' ^^^'^ ^n"' ' respectively. The estimates 
of the S = 1 gap from the QMC method^ and from the DMRG 
one^ are also presented near the left-hand ordinate. 

possible to estimate the gap value irrespective of bound- 
ary conditions. The difference of systematic errors be- 
tween the cases of the periodic boundary condition and 
the twisted boundary condition originates from the char- 
acteristics of the initial sequences. In this meaning, the 
twisted boundary condition is better than the periodic 
boundary condition to estimate the Haldane gap. 

3.2 Cases for S>2 

In this subsection, we estimate the Haldane gap for 
S > 2. We employ the twisted boundary condition to 
obtain the finite-size gap G{N). Since the sequence G{N) 



is increasing with N , it is easy to distinguish whether the 
gap survives or not in the limit of N oo even though 
the gap value is extremely small. 

Let us extrapolate our finite-size gaps for the S = 2 
case by eqs. (6) and (10). The result is summarized in 
Table IV. One finds that data for small N disturb the 
tables In Table IV, A^^'' and data originating from it are 
underlined. If we exclude these underlined data from the 
table, all other data do not disturb the table, which sug- 
gests that the convergence acceleration of the sequence 
starting from iV = 6 is successful. One can also find 
that without the underline show the dependence 

opposite to A^J^\ After excluding the underlined data, 
we do not have a suflacient number of data in B)^' to 

(2) 

know its dependence; thus, we do not employ Bj^g as an 
upper bound in this work. In addition, we cannot con- 
firm whether A^^ and A^g^ are successfully accelerated 
or not, because available ^j^^ originates from Af\ From 
the above reason, it is a careful and reliable judgment 
to consider that A{S = 2) is between A^g^ and B^^ . 
Therefore our conclusion of the estimates for A{S = 2) 
is 

0.0886 ± 0.0018, (12) 

for S = 2. This estimate agrees with A = 0.08917 ± 
0.00004 from the QMC simulation in ref. 6 and A = 
0.0876 ± 0.0013 from the DMRG calculation in ref. 12. 

Next, we study the cases for S > 3. Our numerical re- 
sults under the twisted boundary condition are presented 
in Table V. One can observe that the decay length in 
S* = 4 and 5 is not monotonically increasing. Although 
the decay length in iS = 3 is monotonically increas- 
ing, is negative. This indicates that A^^ and its pos- 
terity should be excluded from the table of acceleration 
according to the above criteria used in the case of S = 2. 
After the exclusion, one cannot judge whether A^^^ is suc- 
cessfully accelerated within the present data. Thus, we 
do not apply the above acceleration procedure to these 
finite-size gaps of 5 > 3 in this study. 

For S > 3, a. serious behavior of the finite-size devi- 
ations appears when one draws a plot of G{N) versus 
the inverse of the system size as a way that is usually 
applied. In this plot, the l/7V-dependence of G{N) re- 
veals concave upwards for small sizes. For S = 3 and 
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Table IV. Sequence of finite-size gaps for the 5 = 2 case under the twisted boundary condition, the convergence acceleration based on 

(k) 

Wynn's algorithm, and upper-bound sequence Bj^ . 



N 




^(0) 


10^ Ay 




-1 ri9 T-i( 1 ) 




iV 


1 r,2 7->(2) 


4 


3.60543243815 
















6 

Q 
O 


5.88636574630 


Z. { Z 














9 










9.046111 








10 


7.57841496.301 


3.30 


8.2962537 












11 










9.125346 








12 


7.93926067768 


3.97 


8.4896523 


4.38 




8.6729401 






13 










9.068141 






8.814100 


14 


8.17386831221 


4.65 


8.6098758 


4.21 




8.7252465 






15 










9.021915 






8.963622 


16 


8.33497991928 


5.32 


8.6881548 


4.66 




8.7779957 


-237.32 





Table V. Finite-size gaps and decay lengths for the 5 > 3 cases under the twisted boundary condition. For 5 = 3, the convergence 
acceleration based on Wynn's algorithm is also applied. 







5 


= 3 




5 = 4 


5 = 5 


N 






lO^A^) 


.(1) 


lO^A^^ 


AO) 


105 A^^ 


AO) 


4 


2.44233786473 








1.29166071630 




0.59777444079 




6 


5.20317763910 








3.52085939468 




2.06995173926 




8 


6.75739139783 


3.48 


8.75932105666 




4.94293026662 


4.45 


3.12913003881 


6.07 


10 


7.66279200087 


3.70 


8.92625320663 




5.80259274045 


3.97 


3.79423190111 


4.30 


12 


8.23027668731 


4.28 


9.18329148633 


-4.63 


6.34716007738 


4.38 







4, the dependence becomes convex upwards for larger 
sizes. For S' = 5, the dependence is stih concave upwards 
in the range up to = 10. If the dependence is concave 
upwards, it is difficult for us to capture a converging be- 
havior. 

Under these circumstances, we take a strategy com- 
posed of the following two steps. The first step is to draw 
a plot of G{N) so that a shape which is concave upwards 
does not appear. The second step is to create a decreas- 
ing sequence from G{N) that is an increasing sequence 
in the new plot. 

In order to carry out the first step, we here introduce 
a renormalized system size N defined as N + Nq so that 
three data for = 4, 6, and 8 reveal a linear dependence 
in the plot of G{N) versus 1/N. The results are depicted 
in Fig. 4. One can observe the 1 / TV-dependence of G{N) 
is always convex upwards in every 5* of Fig. 4 (a)-(d). 
Note in Fig. 4 (a) and (b) that G{N) approaches an 
estimate obtained from QMC in ref. (6) for S = 2 and 3. 
On the other hand, we do not have other estimates with 
which we can compare our G{N) for 5* = 4 and 5. All of 
G{N) in Fig. 4 are increasing with N. Thus, our finite- 
size gaps G{N) under the twisted boundary condition 
are appropriate to estimate the Haldane gaps even for 
S > 2. One finds that G{N) gives lower bounds for the 
Haldane gaps. 

Next, we perform the second step. We focus on neigh- 
boring three data points of system sizes A''— 2, N, and 
iV-l-2 in each panel of Fig. 4. When we apply the fitting 
curve of 



y = C + Dx'' 



(13) 



to the three data points, it is possible to determine the 
parameters C, D, and E uniquely for a given N. Then 
we use G{N), D{N), and E{N) hereafter. Due to the 



above first step, E{N = 6) is necessarily the unity. Note 

that D{N) is monotonically decreasing with increasing 
and that E{N) is monotonically increasing. The result of 
G{N) are also depicted at the corresponding N in Fig. 4. 
One can observe that C{N) is monotonically decreasing 
with increasing A^ in Fig. 4 (a)-(d). At least for 5* = 2 
and 3, G{N) seems to converge from the upper side to 
the gap value estimated from QMC calculations.^ It is 
reasonable to consider that the sequence G{N) becomes 
an upper bound of the gap value. We choose G{N) for 
the largest system sizes as the best upper bound in the 
analysis of Fig. 4 (a)-(d). From the above argument of 
the lower and upper bounds of the gap value, we obtain 

(14) 

(15) 

(16) 



0.0868 ± 0.0034, 



A{S = 2) 
A{S = 3) = 0.0092 ± 0.0010, 
A{S = 4) = 0.00072 ± 0.00009, 



A(5 = 5) = 0.000047 ±0.000010. (17) 

The estimate (14) for 5 = 2 is consistent with the above 
result (12). We can also compare our estimate of A{S = 
3); our estimate agrees with A = 0.01002 ± 0.00003 from 
the QMC simulation in ref. 6. There is no other numeri- 
cal estimate for S* > 4 to the best of our knowledge. Our 
estimate (16) and (17) will be inspected in future if other 
methods become available. 

Note here that the analysis without convergence accel- 
eration is also applicable to the case of 5 = 1; the result 
is 



A{S = 1) = 0.41028 ± 0.0042. 



(18) 



This estimate is consistent with the gap (11). One of the 

differences is that the error of (18) is wider than the one 
of (11). The same situation appears between (14) and 
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Fig. 4. Finite-size gaps G{N) for S = 2, 3, 4, and 5 under the twisted boundary eondition are shown by open cireles. We determine a 
renormahzed system size A'^ defined as Af + A^o so that three data for A' = 4, 6, and 8 reveal a hnear dependence; A^o = —0.3091268, 
1.15226211, 3.047-56114, and 6.25840221 for S = 2, 3, 4, and 5, respectively For S = 2 and 3, the results from QMC in ref. (6) are 
shown by closed squares. Diamonds represent upper bounds C{N). Equation (13) gives C{N) from the finite-size gaps of system sizes 
Ar-2, AT, and N+2. 



(12) in the case of S = 2. These; agreements for 5* = 1 
and 2 suggest that both the methods successfully lead 
to the common estimate that is irrespective of analyzing 
methods. Another difference is that the central value of 
the estimate (18) is slightly smaller than that of the es- 
timate (11). One can also in the case of 5 = 2 observe 
that the central value of the estimate (14) is smaller than 
that of the estimate (12). The estimate from the QMC 
simulation in ref. 6 is larger than the central value of the 
estimate (15) in the case of S' = 3. These facts suggest 
that the present method using eq. (13) may give us the 
central value that is closer to the data of the original se- 
quence. One of the reason may be that C{N) of eq. (13) 
is closer to the true quantity in the thermodynamic limit 
because C{N) is a kind of extrapolated results. It is not 
so easy to know where within the error the true infinite- 
size quantity is only from the data of this work without 
referring other information. However, this question is be- 
yond the present study because a primary purpose of this 
work is the development of a method that makes us pos- 
sible to obtain a reliable systematic error within which 
there exists the true infinite-size quantity. 

Here, we examine the asymptotic formula of the Hal- 
dane gap. 



A(5)=/3|S|2exp(-7r|S|), 



(19) 



for S — > oo from our estimate of the gaps for finite S. 
In order to do it, we introduce new parameters x = 




Fig. 5. Analysis of our estimates of Haldane gaps up to 5 = 5 to 
confirm the asymptotic formula for large S. Errors for 5 = 1 and 
2 are smaller than the symbol size. The dotted line is the one 
with /? ~ 12.8 obtained from the linear fitting of 5 = 4 and 5. 



and y — S~^\og{S'^ / A{S)) when we take the amplitude 
of each spin to be |S| = S in the present analysis. The 
asymptotic formula (19) is rewritten as 



y = IT — X log /?. 



(20) 



Let us input our estimates of the Haldane gaps to A (5) in 
y and plot the x dependence of y. The result is depicted 
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in Fig. 5. One can find a linear behavior for finite but 
large S up to S = 5. Wc have fitted our data for S = 4 
and 5 with the straight line (20); the best fit is produced 
by /3 = 12.8 ± 1.5. The linear behavior suggests that the 
asymptotic formula (19) holds well for large S. Note here 
that if one uses the formula (19) and our estimate of (3, 
one can predict the magnitude of the Haldane gap for 
,S > 5. 

4. Summary and Remarks 

We have developed methods to estimate a systematic 
error in the extrapolation of finite-size data obtained 
from numerical diagonalizations of small clusters. The 
methods have been applied to study the Haldane gaps 
of the one-dimensional integer-S' Heisenberg antiferro- 
magnet. We have demonstrated that the methods work 
well to the finite-size data under the twisted bound- 
ary condition. Our best estimate of the 5=1 case is 
A = 0.4104789 ±0.0000013. We are also successful in ob- 
taining the gaps up to 5 = 5 even though the magnitudes 
of the gaps are extremely small. Our estimates of the gaps 
for large S make us confirm that the asymptotic formula 
of the gap holds well. We have found that the twisted 
boundary condition is more appropriate than the peri- 
odic and open boundary conditions in order to extrapo- 
late finite-size data of the Haldane gaps to the thermo- 
dynamic limit. However, the most appropriate boundary 
condition is not always the twisted boundary condition. 
The most appropriate boundary condition depends on 
the model and physical quantities. The present work 
strongly suggests that the examination of various bound- 
ary conditions is useful to know reliable quantities in the 
thermodynamic limit. It should be examined in future 
studies whether the methods work or not in extrapola- 
tions of other quantities, like the correlation functions 
of a system. The present method could give detailed in- 
formation about the quantum state, which contributes 
much to a deeper understanding of properties of the sys- 
tem. 

Finally, it is noticeable that to produce a new sequence 
by eq.(lO) is a quite versatile way. Only Wynn's transfor- 
mation is employed in this work, but eq.(lO) is available 
irrespective of acceleration methods. The parameter of 
sequences is not limited to the system size. Usefulness 
of this way for other parameters should be examined. In 
principle, the way of B}^' is applicable when the num- 
ber of data in an original sequence A^'' is five at least 

because the direction of the dependence of B^j^^ can be 
checked. Note also that properties of the new sequence 
B^^ are related to whether the acceleration procedure is 
successful or not as we have mentioned in Appendix. 
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Appendix: Convergence acceleration and new 
sequence 

Let us consider monotonic scalar sequences T*^ and Sn 
that share the common limit S. According to ref. 17, 
the convergence of r„ is faster than that of Sn when the 
following condition is satisfied 



lim 



= 0. 



(A-l) 



n— >oo Sn — S 

In a practical problem, however, lengths of the sequences 
are finite; S is not a known quantity but the one that 
should be obtained. Therefore, the condition (A-l) can- 
not be examined directly. A substitute condition for 
(A-l) is given by 



lim 

n— >(X) tin 



■ T 



= 0. 



(A-2) 



for an arbitrary integer m(> 0). eq. (A-2) suggests that 



1 > 



> 



T —T 



Sn-\-5n Sn-\-Sn-\-m 



'->n+25n ~ '-'n+25n+m 

for n in the region where a converging behavior appears 
in {Tn—Tn+m)/{Sn—Sn+m)- Notc here that this equation 
with dn = m = 2 is related to eq. (9). 

Let us get back to eq. (A-l). Suppose that both Sn and 
T„ approach S from the same side. The equation (A-l) 
suggests that 

^"+'" ~ ^ < (A-4) 

i~>n+m ~ i~> i~>n ~ iJ 

for sufficiently large n and positive m. This inequality is 
rewritten as 

[{Sn+m Sn^ (^n+m '^n)]S <^ TnSn+vn '^n+mSm 

(A-5) 

because Sn is monotonic. When T„ is obtained from Sn 
through a successful acceleration procedure and when Sn 
and Tn are monotonically increasing, one can find 

'^nSn+m T-^n+mSn 



s < 



(A-6) 



The right hand side of this inequality can be written to 
be 



S + 



Tn-S 



S„4 



{Sn - S) 



1 



(A-7) 



which is easily found to converge to S in the limit of 
n — > DO with a help of the limit (A-2). 



Let us consider the case of r„ 



A. 



(fe) 



Sn 



A 



(k-i) 



and 



: 2: The right hand side of the inequality (A-6) 
becomes b':^]^^. From the inequality (A-6) and the limit 

of (A-7), b'^^1^ defined as eq. (10) is found to converge 
to S from the upper side. Note that the direction of the 
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Table A-1. An application of the convergence acceleration to the sequence (A-8) having its limit as tt. The decay lengths obtained from 



each sequence by cq. (7) arc also acompanied. 
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i V 


AO) 


AO) 
?/v 




.(1) 


4(2) 




4(3) 


A3) 




.(4) 


4(5) 
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A 
























Q 

8 


3.106348882832 


2.48 
2.97 


3.12083429 
3.13238707 


















10 


3.120142116779 


3.37 


3.13715611 


2.26 


3.13939334 














12 


3.128166065374 


3.69 


3.13932581 


2.54 


3.14068137 














14 


3.133009091490 


3.96 


3.14038271 


2.78 


3.14118662 


2.14 


3.14135408 










16 


3.136013726258 


4.19 


3.14092445 


2.99 


3.14140146 


2.34 


3.14149788 










18 


3.137917979434 


4.39 


3.14121336 


3.18 


3.14149861 


2.52 


3.14155251 


2.07 


3.14156654 






20 


3.139145542642 


4.56 


3.14137243 


3.35 


3.14154473 


2.68 


3.14157476 


2.23 


3.14158257 






22 


3.139947946854 


4.70 


3.14146233 


3.50 


3.14156750 


2.83 


3.14158434 


2.37 


3.14158853 


2.02 


3.14158978 



(k) 

convergence of B^_(_-^ is opposite to that of T„ and that 
of 5„. 

The above argument is appUcable when 5„ and T„ are 

(k) 

monotonically decreasing. In this case, -B„_|_i converges to 
S from the lower side. The direction of the convergence of 
b'^^-^ is opposite to that of T„ and that of Sn irrespective 
of the direction of T„ and Sn- 

Finally, a numerical example is presented for a demon- 
stration of the procedure proposed in this paper. We con- 
sider the sequence {A^^^} = {A^°\ Af\ a'^\ A^°\ ■ ■ ■ } 
defined as 



.(0) 



fe=o 



{2k)\ 



2fc + r 2 ^ 



(A.; 



where {2k)U = (2fc)(2fc - 2) • • • 4 • 2 and {2k - 1)!! = 
{2k — l)(2fc — 3) • • • 3 • 1. This initial sequence is obvious 
to increase monotonically. It is easily understood that 
the limit of TV ^ 00 is given by 3 arcsin(\/3/2) = 7r(= 
3.14159265- • •) if one remembers the Taylor exapansion 
of arcsina; in a neighborhood of a; = 0. If we apply 
Wynn's transformation (6) and create a new sequence 



B 



(fe) 



N+ 



aP and 

. The result is summarized in Table A-1, 



J by means of eq. (10) in the case of Tjv 

' (k-l) 



Table A- 2, and Fig. A-1. One can confirm a successful 

(k) 

acceleration of convergence of A)^ in Table A-1 from a 
judgement based on Conditions I, II, and III in §3.1.2. 
A successful creation of B^^ which decreases monoton- 
ically with respect to TV can be observed in Table A-2. 
Figure A-1 is presented in order to confirm that the fea- 
tures of A^ and hold for even larger N. One can 



find that each approaches the exact value from the 

t (k) 

upper side while A)^ approaches it from the opposite, 
namely lower, side. It is reasonable to consider that each 



datum of A^ and B)!^' is a lower bound and a upper 
bound, respectively, for the rigorous limit n. 
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Fig. A-1. Log-log plots of (a) n — A 



(fe) 



and (b) 



n of se- 



quence (A-8) for N up to 30. Reversed triangles, circles, squares. 



triangles, and crosses in (a) denote ir — , tt — A 



l(2) 



and tt- 



(4) 



respectively. Circles, squares, triangles. 



and crosses in (b) denote — tt, B 



,(2) 



TT, — TT, and 



TT, respectively. 
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